Anisotropic excitation spectrum of a dipolar quantum Bose gas 
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We measure the excitation spectrum of a dipolar Chromium Bose Einstein Condensate with 
Raman-Bragg spectroscopy. The energy spectrum depends on the orientation of the dipoles with 
respect to the excitation momentum, demonstrating an anisotropy which originates from the dipole- 
dipole interactions between the atoms. We compare our results with the Bogoliubov theory based 
on the local density approximation, and, at large excitation wavelengths, with numerical simulations 
of the time dependent GrossPitaevskii equation. Our results show an anisotropy of the speed of 
sound. 
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Interactions play a major role in the physics of Bose 
Einstein Condensates (BECs) made of trapped neutral 
atoms. Attractive interactions may lead to collapse 
while repulsive interactions confer a collective nature to 
excitations Q, and lead to superfluidity Q. In the first 
produced BECs only contact interactions played a signif- 
icant role. The study of long-range anisotropic (partially 
attractive) dipole dipole interactions (DDIs) in quantum 
gases was opened by the production of Chromium BECs 
y , followed more recently by experiments with Dyspro- 
sium [s'j'i^ and Erbium Q. DDIs introduce an anisotropy 
in the expansion dynamics Q and the collective exci- 
tations of a Cr BEC Q. Moreover when DDIs over- 
whelm contact interactions a BEC undergoes a charac- 
teristic d-wave like implosion which reveals the structure 
of DDIs 0, [lol |. In addition, anisotropic dipolar interac- 
tions may lead to the study of anisotropic superfluidity in 
Bose 11] or Fermi gases [1^, with possible analogies to 
anisotropic superconductivity in cuprates. In this paper 
we study elementary excitations of a Cr BEC and show 
an anisotropic excitation spectrum, providing a signa- 
ture of an anisotropy of the speed of sound and there- 
fore, according to the Landau criterion, a possibility for 
anisotropic superfluidity. 

For investigating elementary excitations with a well de- 
flned momentum, Raman-Bragg spectroscopy has devel- 
oped into a powerful instrument [l3| . Following the first 
seminal demonstration that rapidly followed the creation 
of atomic BECs [3, 0, more systematic series of ex- 
periments gave a full picture of the BEC excitation spec- 
trum from the low-energy phonon-like excitations to the 
high-energy single-particle regime 16[. Despite several 



reviews and theoretical papers pointing out the inter- 
est of the excitation spectrum of bosonic gases featuring 
DDIs 17-l3|, experimental results are not yet available 
(ioj . Chromium atoms away from a Feshbach resonance 
are particularly suitable to study bulk properties of dipo- 
lar BECs since DDIs are non negligible without leading 
to collapse. In this work, we provide the first evidence 
that long-range dipolar interactions induce an anisotropy 



of a BEC excitation spectrum using Raman-Bragg spec- 
troscopy. The spectrum is anisotropic as the resonance 
condition depends on the angle, 9, between the polariza- 
tion axis, parallel to the magnetic field B, and the excita- 
tion wave- vector q (see Fig[T]). We probe this anisotropy 
throughout the whole excitation spectrum with differen- 
tial measurements for two orthogonal orientations of the 
magnetic field, either parallel or orthogonal to q. 
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FIG. 1: (Color online) Principle of the experiment and absorption 
pictures after a Bragg pulse. The two Bragg beams are coupled to the 
BEC with an angle 0, and transfer a momentum ?iq and an energy hf 
to the excited fraction. Due to DDIs, the excitation spectrum depends 
on the angle 9 between q and B. The two absorption pictures show in 
false color the density after Bragg pulse and time of flight, and solid 
lines indicate the direction of momentum transfer. In a) the small value 
of (p (14deg), and hence of g, barely allows spatial separation of the 
excited fraction, contrary to the case of large (p (82 deg) in b). 



The anisotropy of the excitation spectrum for an ho- 
mogenous dipolar BEC is well understood within the Bo- 
goliubov theory (l7l |: 
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where e(q) and q are respectively the energy and the 
wave vector of the excitation, m the atom mass, n the 
atomic density, and g = Airh^a/m (with a the scat- 
tering length); €dd is the dimensionless parameter scal- 
ing the relative importance of dipolar interactions with 
respect to contact interactions. For magnetic DDIs, 
^dd = /io/^m™/12?i^a (/iQ is the vacuum permeability, /z^ 
the atom magnetic moment equal to 6 Bohr magnetons 
for Cr). In absence of DDIs {edd = 0), one recovers the 
known spectrum for a BEC with only contact interactions 
m. For Cr edd is 0.16 (with a = 102. 5ao p3|), which al- 
lows to explore the anisotropic character of eq. ^ more 
easily than for alkali BECs (e.g. edd = 0.01 for Rb). 
The angular dependance of eq. ([T]) can be interpreted 
as an attractive contribution (resp. repulsive) of DDIs 
to the excitation energy for the perpendicular (parallel) 
case p^. It is a direct consequence of the momentum 
sensitivity of DDIs, in contrast to contact interactions 
which are momentum-insensitive. 

According to eq.([T]), the sound velocity defined as 
cg = limq_j.o^~^; becomes anisotropic. It is maximal 
in the parallel {9 — 0) geometry, while it is minimal 
in the perpendicular {9 = 7r/2) geometry. Defining 
the velocity without DDIs cq — {gn/mY/'^ ^ one obtains 
C|| = co(l -I- 2edd)^/^, and c_l = co(l - Cdd)^^^, so that 
C||/c_L = 1-25 for Cr. This dipolar effect is larger than 
previous observations [§, 0], where effects are in the few 
% range, much less than edd, as a result of an angular 
averaging of DDIs in the BEC [1, 23|. The situation is 
quite different for the BEC excitations described by eq. 
([J), as there is no such angular averaging: 6 has the same 
value for all the atoms. 

To obtain the experimental excitation spectra, we first 
create a chromium BEC in a crossed-beam dipole trap. 
Details on the experimental setup and the procedures 
can be found in 2J]. The BEC, comprising about 10000 
atoms polarized in the absolute ground (Zeeman) state, 
is confined at the bottom of the trap with frequencies 



y, 



= 27r X (145, 180, 260) Hz, within a magnetic 
field of B=50 mGauss. To impart a momentum g to a 
fraction of the condensed atoms, we use two intersecting 
phase-locked focused laser beams derived from the same 
single- mode solid-state laser (A = 532 nm) . Two acousto- 
optic modulators (AOMs) shift the beam frequencies by 
respectively and D, + 2nf, with = 27r x 80 MHz. At 
the atom location, their respective waists are 40 and 32 
/im while their powers are between 500 and 1000 iiW. The 
beam intensities are chosen to limit the excited fraction 
to about 15 % which ensures a good trade-off between 
the validity of the perturbative Bogoliubov approach for 
theoretical interpretation and the signal to noise ratio. 
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FIG. 2: (Color online) Excitation spectra for g^o — 0.8 (0 — 14deg). 
The excited fraction Pcxc is plotted versus the detuning frequency / be- 
tween the two Bragg beams, for ^ — (filled circles) and 9 — 7r/2 (open 
circles). The two solid lines are asynietric gaussian fits to the data (see 
text). The error bars represent the 1 sigma statistical uncertainty as- 
sociated to three measurements. The two spectra are shown separately 
below ( a): 9 — 7r/2, b): 9 — 0) along with the results of our linear 
response theory (dashed lines) which has no adjustable parameters. 



The angle (j) between the beams propagating directions 
sets the momentum value q {q = 2'k/Xx sin((/)/2)), while 
the energy difference hf, equal to the energy of the ex- 
citation, is adjustable at will. Although optical access 
restricts the accessible values for (p, thus also that for q, 
we have been able to probe the dispersion relation from 
the phonon to the free particle regime. 

The two beams hereafter referred to as the Raman 
beams are switched-on only for a brief pulse-time. Ra- 
man transitions are efficiently driven between the rest 
ground state and excited states with momentum q and 
energy hf when the beam photon energy difference 
matches the excitation energy i.e. when the Raman- 
Stokes resonance condition is met. Since the process can 
be interpreted as diffraction of the matter-wave onto the 
moving lattice created by the Raman-beams, the whole 
process is often referred to as Bragg spectroscopy of the 
excitation spectrum of the BEC. To ensure that the mo- 
mentum of the atoms is solely set by the Bragg pulse 
(and is not modified by trap dynamics), the pulse dura- 
tion TBragg should be far less than the oscillation periods 
in the trap IJ] . Meeting this requirement causes a strong 
Fourier broadening and we found a trade-off by setting 



3 



TBragg = 1-5 ms: the Fourier broadening (half width at 
l/e^=300 Hz) is then smaUer that the experimental spec- 
tra widths, and spatial separations between the excited 
fraction and the ground state are close to the ones ex- 
pected for the Bragg momentum transfer (see Fig[T]). 

We measure the atom momentum distribution by re- 
leasing the atoms from the trap after the Bragg pulse 
and by performing absorption imaging after an expan- 
sion time of 5 ms (see Fig [ij . From the density profiles 
of these distributions along the excitation direction, we 
infer the excited fraction Pcxc = [-^cxc/(-^exc+-^o)] where 
A'cxc is the number of excited atoms and Nq the number 
of atoms remaining in the ground g ~ state. We record 
Bragg spectra by monitoring the excited fraction versus / 
for a given value of (p. It is clear from Fig [2] that the exci- 
tation energy spectrum for parallel polarization is shifted 
towards high frequency with respect to the spectrum for 
orthogonal polarization, as expected from eq. ([T]). 

The results of eq. ([1]) hold only in the homogeneous 
case. With trapped gases, two kinds of finite size effects 
have to be taken into account, which are sources of spec- 
tral broadenings: the density inhomogeneity in the trap, 
and the non zero-width momentum distribution (induc- 
ing Doppler effect), which are dominant at respectively 
low and high excitation momentum [25'|. To account for 
our experimental results, we extended the theory devel- 
oped in [2^ based on the local density approximation 
(LDA), by including DDIs. The validity of LDA with 
DDIs, discussed in [18|], is guaranteed if g 3> l/RrFmin 
with RrFmin the smallest BEC Thomas Fermi (TF) 
radius. Then the BEC can be considered as a locally 
homogeneous 3D system [27[. q 3> i/RrFmin also en- 
sures that effects of energy discretization can be safely 
neglected. 

To compare with our experimental spectra (see Fig 
[21), we take into account the Doppler width and make 
a convolution between the corresponding Gaussian and 
the LDA excitation spectrum [28 1. The Fourier broad- 
ening due to the finite excitation duration is also taken 
into account within the linear response theory (29| . The 
two photon Raman frequency which sets the excitation 
amplitude |16l | is estimated from a calibration of the lat- 
tice depth created by the Raman beams . The agree- 
ment with theory (with no adjustable parameter) for the 
spectra line shapes is good except for the two lowest ex- 
perimental values of q. For the excitation amplitude, the 
discrepancy remains at most in the 50% range for all 
spectra. All features are in very good agreement in the 
case of Fig O 

We have performed a study of effects of DDIs by prob- 
ing the excitation spectrum as a function of q from the 
phonon regime to the free particle regime (see Fig[3]). We 
define the experimental value of e(q) as the mean value 
of the excitation spectra, obtained by using asymmetric 
Gaussian fits; this procedure empirically takes into ac- 
count the non-symmetric nature of the excitation spec- 
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FIG. 3: (Color online) a) Energy spectra: the energy €(q) divided 
by h is plotted as a function of Bragg-transfcrrcd momentum q. Cir- 
cles correspond to the experimental data: for a given value of g, the 
top (resp. bottom) circle corresponds to the parallel (perpendicular) 
case. The point at q^o — 1.6 corresponds to a four photon process (see 
text). The two solid lines are the results of LDA calculations for resp. 
parallel (top line) and orthogonal (bottom line) case. The dashed line 
gives the results for a free particle. Inset: zoom for low q values, b) 
Relative shift A{q) of the excitation energies (see text) as a function 
of q. Black points: results of fits of the experimental data, with er- 
ror bars. Solid line: results of LDA calculations. The diamonds show 
results of numerical simulations. For a) and b), vertical error bars rep- 
resent IfT statistical uncertainty, while horizontal ones correspond to 
uncertainties on values of (p. 



tra which is expected from the Thomas Fermi distribu- 
tion (see for example ^^). We define e\\{q) (respectively 
e±{q)) the mean value obtained for the parallel (orthogo- 
nal) case. Fig [3^) shows the corresponding data. We use 
the DDI-independent healing length — 7i/(2m^no)^/^ 
(with no the measured BEC peak density) to normal- 
ize the q axis: indicates the frontier between the 
phonon domain (g^o < 1) ^^^d the single-particle domain 
{q^o ^ !)• The data correspond to four (j) angles. To 
fill the gap between the accessible low (6°, 10° and 14°) 
and the high (82°) values of 0, we realized a four-photon 
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Raman-like excitation [3l| to obtain data for another q 
value {q = 47r/A x sin(1472), and q^o = 1.6). The hori- 
zontal error bars are mainly set by the uncertainty on the 
experimental value of and the vertical ones are deter- 



mined by applying a bootstrapping procedure [32[ to the 
experimental data. We compare predictions of our the- 
ory assuming a constant (10^) number of atoms with the 
corresponding experimental data in Fig [3^) [11] . We ob- 
tain a good general agreement between data and theory 
with no free parameter; some of the discrepancies may 
originate from an underestimation of the uncertainties 
on the values of (p. 

To emphasize how the two dispersion curves differ, we 
plot in Fig [Sb) the relative variation of e(q), A(g) = 
measure an anisotropy of the excita- 
tion spectrum of the Cr-BEC, both in the low q limit, 
and the high q limit. As systematic effects on /S.{q) re- 
lated to changes in trapping frequencies and atom num- 
ber induced by the change in 9 remain negligible (from 
at most 1.5% at low q to below 0.1% at high g), our 
results demonstrate a DDI induced anisotropy of the ex- 
citation spectrum. The agreement with our linear re- 
sponse theory is relatively good, except for the lowest 
value of q. This is not surprising since then the excita- 
tion wavelength is larger than the BEC Thomas Fermi 
diameter (equal to 7 fim). For long wavelengths, discrete 
modes can be excited, for which the effects of DDIs are 
much smaller than for phonons, as observed in Q. We 
have therefore solved numerically the time dependent 3D 
Gross-Pitaevskii equation including contact interactions 
and DDIs. Our numerical simulations are in rather good 
agreement with our experimental data (see FiglSb)), and 
show that relative-shifts significantly smaller than the 
ones obtained with LDA theory are expected as the ex- 
citation wavelength increases. 

We finally rely on LDA theory to extract quantita- 
tive values for the sound velocity from experimental data. 
In our case the extensive experimental study of the lin- 
ear phonon part of the spectrum is not accessible due 
to the small size of our BEC (contrary to [l3|)- We 
therefore chose to derive sound velocities from the data 
at 9^0 = 0.8, as this is the lowest value of q for which 
we find a good agreement with LDA theory. The corre- 
sponding excitation frequency (around 1.5 kHz) is signif- 
icantly higher than the trap frequencies, which explains 
why energy discretization effects are small. From the ex- 
perimental values e\\/h = (1.57 ± 0.05) kHz and e^/h ~ 
(1.35 ± 0.03) kHz, we derive sound velocities through the 
linear response theory. The Feynman law (see jig ]) re- 
lates the energy gq to the (LDA averaged) static struc- 



ture factor 5a 



with 



h q /2m the free 



particle energy; and is related to eg through 

15 / 3+aq _ (3+2aq-aq^) ' 
4 I 



mm.s ^ and cy = 2.64 ±0.1 mm.s ^ , in good agreement 
with LDA predictions in our experimental conditions 
equal to 2.02 ± 0.05 mm.s"^ and 2.53 ± 0.05 mm.s'^ re- 
spectively {cg^LDA = iIf [^(1 + erf<i(3cos2 - 1))]'/'). 

In this article, we have shown that the excitation spec- 
trum of a spin-polarized chromium BEC is anisotropic as 
a consequence of the dipolar character of the interactions 
between the atoms. As long as the wave vector is not too 
small compared to the inverse of the BEC size, our ex- 
perimental results are correctly accounted for by a linear 
response theory based on the Bogoliubov approach. The 
demonstrated existence of an anisotropic speed of sound 
raises the question of anisotropic superfluidity, for exam- 
ple how excitations in different directions couple. Exten- 
sion of our work to BECs with larger DDIs could allow to 
discover rotonic features in the excitation spectrum 
of quasi 2D BECs. 
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We obtain c_l = 2.06 ± 0.05 
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